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Abstract 

The angular distribution and differential branching fraction of the decay 
2?^— 7- K^fi^n'^ are studied with a dataset corresponding to 1.0 fb^^ of integrated 
luminosity, collected by the LHCb experiment. The angular distribution is mea- 
sured in bins of dimuon invariant mass squared and found to be consistent with 
Standard Model expectations. Integrating the differential branching fraction 
over the full dimuon invariant mass range yields a total branching fraction of 
B{B+^ K+fi+fj,-) = (4.36 ±0.15 ±0.18) x lO"'^. These measurements are the most 
precise to date of the — K~^fj.'^^~ decay. 
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1 Introduction 



The — jj,^ jj,^ deca}|j is a 6 — > s flavour changing neutral current process that is 
mediated in the Standard Model (SM) by electroweak box and penguin diagrams. In 
many well motivated extensions to the SM, new particles can enter in competing loop 
diagrams, modifying the branching fraction of the decay or the angular distribution of the 
dimuon system. The differential decay rate of the (B^) decay, as a function of cos6i, 
the cosine of the angle between the fi~ (/i"^) and the {K~) in the rest frame of the 
dimuon system, can be written as 

I dV\B+ ^ K+ a+ a-] 3, n,,l 

which depends on two parameters, the forward-backward asymmetry of the dimuon system, 
A-pB-i and the parameter Fh (i]|2). If muons were massless, Fh would be proportional to the 
contributions from (pseudo-) scalar and tensor operators to the partial width, F. The partial 
width, Afb and Fh are functions of the dimuon invariant mass squared (g^ = 

In contrast to the case of the B^ — )■ K*^fi~^n~ js] decay, Aps is zero for B^ — )■ K^jj,^jj,~ 
in the SM. Any non-zero value for ApB would point to a contribution from new particles 
that would extend the set of SM operators. In models with (pseudo-)scalar or tensor-like 
couplings \Apb\ can be enhanced by up to 15% |2]|4|. Similarly, Fu is close to zero in 
the SM (see Fig. |3]), but can be enhanced in new physics models, with (pseudo-) scalar or 
tensor-hke couplings, up to Fh < 0.5. Recent predictions for these parameters in the SM 
are described in Refs. (2||5}[6]. Any physics model has to satisfy the constraint |Afb| < -Ph/2 
for Eq. ([T]) to stay positive in all regions of phase space. The contributions of scalar and 
pseudoscalar operators to ApB and Fh are constrained by recent limits on the branching 
fraction of B^ — )■ [7] . The differential branching fraction of B~^ — )■ ^~ can 

be used to constrain the contributions from (axial-) vector couplings in the SM operator 
basis |6)|8)|9). 

The relative decay rate of B^ — )■ K~^^'^fi~ to B^ — )■ K'^^~^fi~ has previously been studied 
by the LHCb collaboration in the context of a measurement of the isospin asymmetry |10| . 
This paper presents a measurement of the differential branching fraction (di3/dg^), Fh 
and ApB of the decay B~^ — )■ K~^ii~^^~ in seven bins of and a measurement of the total 
branching fraction. The analysis is based on 1.0 fb~^ of integrated luminosity collected in 
^/s = 7TeV pp collisions by the LHCb experiment in 2011. 

2 Experimental setup 

The LHCb detector |lli| is a single-arm forward spectrometer, covering the pseudorapidity 
range 2 < r/ < 5, that is designed to study b and c hadron decays. A dipole magnet with a 
bending power of 4 Tm and a large area tracking detector provide a momentum resolution 
ranging from 0.4% for tracks with a momentum of 5 GeV/c to 0.6% for a momentum of 

^Charge conjugation is implied throughout this paper unless explicitly stated otherwise. 
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100 GeV/c. A silicon micro-strip detector, located around the pp interation region, provides 
excellent separation of B meson decay vertices from the primary pp interaction and an 
impact parameter resolution of 20 |xm for tracks with high transverse momentum (pt)- Two 
ring-imaging Cherenkov (RICH) detectors provide kaon-pion separation in the momentum 
range 2 — lOOGeV/c. Muons are identified based on hits created in a system of multiwire 
proportional chambers interleaved with iron filters. The LHCb trigger comprises a hardware 
trigger and a two-stage software trigger that performs a full event reconstruction. 

Samples of simulated events are used to estimate the contribution from specific sources 
of exclusive backgrounds and the efficiency to trigger, reconstruct and select the — t- 



K'^^^fi signal. The simulated pp interactions are generated using Pythia 6.4 12 with 



a specific LHCb configuration |13|. Decays of hadronic particles are then described by 



EvtGen |14| in which final state radiation is generated using Photos 15 . Finally, the 
Geant4 toolkit |16| is used to simulate the detector response to the particles produced 
by Pythia/EvtGen, as described in Ref. [17]. The simulated samples are corrected for 
differences between data and simulation in the momentum spectrum, the detector 
impact parameter resolution, particle identification and tracking system performance. 



3 Selection of signal candidates 

The B^ — )■ K^fi^^'^ candidates are selected from events that have been triggered by a 
single high transverse- momentum muon, with > 1.5GeV/c, in the hardware trigger. In 
the first stage of the software trigger, candidates are selected if there is a reconstructed 
track in the event with high impact parameter with respect to the primary pp interaction 



and high px 18 . In the second stage of the software trigger, candidates are triggered on 



the kinematic properties of the partially or fully reconstructed B^ candidate 19 . 

Signal candidates are then selected for further analysis based on the following re- 
quirements: the B^ decay vertex is separated from the primary pp interaction; the B~^ 
candidate impact parameter is small, and the kaon and muon impact parameters large, 
with respect to the primary pp interaction; the B~^ candidate momentum vector points 
along the B~^ line of flight to one of the primary pp interactions in the event. 

A tighter multivariate selection, using a Boosted Decision Tree (BDT) j20] with the 



AdaBoost algorithm 21 , is then applied to select a clean sample of B^ — )■ K^fi'^fi~ 
candidates. The BDT uses kinematic variables including the reconstructed B^ decay 
time, the angle between the B~^ line of flight and the B~^ momentum vector, the quality 
of the vertex fit of the reconstructed B~^ candidate, impact parameter (with respect to 
the primary pp interaction) and pt of the B~^ and muons and the track quality of the 
kaon. The variables that are used in the BDT provide good separating power between 
signal and background, while minimising acceptance effects in and cosdi that could 
bias the differential branching fraction, Apsici'^) or FH(g^). The multivariate selection 
is trained on data, using B~^ — )■ J/ip [Jjip — )■ /i^^/i") candidates as a proxy for the 
signal and B^ K^fi~^fi~ candidates from the upper mass sideband (5350 < mK+^+^- < 
5600MeV/c^) for the background. The training and testing of the BDT is carried out 
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using a data sample corresponding to 0.1 fb^^ of integrated luminosity, that is not used in 
the subsequent analysis. The BDT selection is 85 — 90% (depending on q^) efficient on 
simulated candidates that have passed the earlier selection. 



Finally, a neural network, using information from the RICH 22 , calorimeters and 
muon system is used to reject backgrounds where a pion is incorrectly identified as the 
kaon from the — )■ K~^n~^ij,~ decay. The network is trained on simulated event samples 
to give the posterior probability for charged hadrons to be correctly identified. The particle 
identification performance of the network is calibrated using pions and kaons from the 
decay chain t- Z)°(— )■ K~'k~^)7i~^ in the data. Based on simulation, the efficiency of 

the neural network particle identification requirement is estimated to be > 95% on the 
signal. 

The contribution from combinatorial backgrounds, where the reconstructed K^, /x"*" 
and fi~ do not come from the same 6-hadron decay, is reduced to a small level by the 
multivariate selection. Remaining backgrounds come from exclusive 6-hadron decays. 
The decays -B"*" — ?• K^J/ip and ^ K~^iIj{2S) are rejected by removing the regions of 
dimuon invariant mass around the charmonium resonances (2946 < m^+^- < 3176MeV/c^ 
and 3586 < m^+^- < 3776MeV/c2). Candidates with ttik+^+i,- < 5170MeV/c2 were also 
removed to reject backgrounds from partially reconstructed B decays, such as B^ — t- 
K*^fi^fi''. The potential background from t- K^J/ip [J/ip — )■ fi'^fi~), where the kaon 
is identified as a muon and a muon as the kaon, is reduced by requiring that the kaon 
candidate fails the muon identification criteria if the K~^fi~ mass is consistent with that 
of the J/ip or iIj{2S). Candidates with a K^fi~ mass consistent with coming from a 
misidentified Z}°— t- K^tt^ decay are rejected to remove contributions from 5+— > D^n^. 
After the application of all of the selection criteria, the dominant sources of exclusive 
background are S^— ?• K^tt^-k^ \23\ and B^ ^ 7r+/i+/x~ [24j. These are determined from 
simulation to be at the level of (1.5 ± 0.7)% and (1.2 ± 0.2)% of the signal, respectively. 



4 Differential and total branching fraction 

The K+yU^yU^ invariant mass distribution of the selected B^ — )■ candidates is 

shown in Fig. [1} The number of signal candidates is estimated by performing an extended 
unbinned maximum likelihood fit to the K~^fi'^fi~ invariant mass distribution of the 
selected candidates. The signal line-shape is extracted from a fit to a B~^ — )• K~^J/ijj 
{J/ip — )■ control sample (which is two orders of magnitude larger than the signal 



sample), and is parameterised by the sum of two Crystal Ball functions 25 . The 
combinatorial background is parameterised by a slowly falling exponential distribution. 
Contributions from B^ — t- K^tx^tt~ and B^ — > 7r+/i+yU^ decays are included in the fit. 
The line shapes of these peaking backgrounds are taken from simulated events. In total, 
1232 ±40 B+^ K+yi+^- signal candidates are observed in the 0.05 < < 22.00 GeVVc^ 
range. The yields in each of the (f bins used in the subsequent analysis are shown in 
Tabled 

The differential branching fraction in each of the q^ bins is estimated by normalising the 
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Figure 1: Invariant mass of selected — t- K^fi^fi'^ candidates with 0.05 < < 
22.00 GeV^/c^. Candidates with a dimuon invariant mass consistent with that of the 
J/ip or ip{2S) are excluded. The peaking background contribution from the decays 
B'^ — > fC+vr+vr" and B~^ — )• tt^ ^~ is indicated in the figure. 



" Theory ^■Binned theory 

-•-LHCb 

T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




[GeV^/c^] 



Figure 2: Differential branching fraction of ?■ K^n^n" as a function of the dimuon 
invariant mass squared, g^. The SM theory prediction (see text) is given as the continuous 
cyan (light) band and the rate-average of this prediction across the bin is indicated by 
the purple (dark) region. No SM prediction is included for the regions close to the narrow 
cc resonances. 
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— 7- K^fi^fi event yield, A^sig, in the bin to the total event yield of the — )■ K^J/ip 
sample, N^+j/^, and correcting for the relative efficiency between the two decays in the q"^ 
bin, eK+ji-iii/sK+ti+iJ.--, 



The branching fractions of B^ — )■ K^J/ip and J/?/' — > are B{B^ — )■ K^J/ip) = 

(1.014 ± 0.034) X 10-3 and i3(J/V' = (5.93 ± 0.06) x 10^^ [26] . The resulting 

differential branching fraction is shown in Fig. |2} 

The bands shown in Fig. [2] indicate the theoretical prediction for the differential 
branching fraction and are calculated using input from Refs. |6] and [27] . In the low 
region, the calculations are based on QCD factorisation and soft collinear effective theory 
(SCET) [28], which profit from having a heavy B^ meson and an energetic kaon. In the 
soft-recoil, high region, an operator product expansion (OPE) is used to estimate the 
long-distance contributions from quark loops [29||30]. No theory prediction is included in 
the region close to the narrow cc resonances (the J/ip and il!{2S)) where the assumptions 
from QCD factorisation/ SCET and the OPE break down. The form-factor calculations 
are taken from Ref. j3l]. A dimensional estimate is made on the uncertainty on the decay 
amplitudes from QCD factorisation/SCET [32] . 

The total branching fraction is measured to be 

B{B+^ K-^fx+f^-) = (4.36 ± 0.15 ± 0.18) x 10"^ , 

by summing over the partial branching fractions and accounting for the regions that 
are not used in the differential branching fraction analysis. These regions account for 
~ 14.3% of the total branching fraction (no uncertainty is assigned to this number). This 
estimate ignores long distance effects and uses a model for dF/dg^ described in Ref. [l] to 
extrapolate across the cc resonance region. The values of the Wilson coefficients and the 
form-factors used in this model have been updated according to Refs. 133) and [sT]. 



5 Angular analysis 

In each bin of g^, Ape and are estimated by performing a simultaneous unbinned 
maximum likelihood fit to the K^fx^fi'^ invariant mass and cos^£ distribution of the 
B~^ candidates. The candidates are weighted to account for the effects of the detector 
reconstruction, trigger and the event selection. The weights are derived from a SM 
simulation of the -B"*" — K~^fx~^fi~ decay in bins of width 0.5GeV^/c^ in g^ and 0.1 in 
cosOe- This binning is investigated as a potential source of systematic uncertainty. The 
largest weights (and largest acceptance effects) apply to events with extreme values of 
cos6'£ (I cos 6*^1 ~ 1) at low g^. This distortion arises mainly from the requirement for a 
muon to have p > 3 GeV/ c to reach the LHCb muon system. This effect is well modelled 
in the simulation. 
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Figure 3: Dimuon forward-backward asymmetry, Ap^, and the parameter for t- 
function of the dimuon invariant mass squared, q^. The SM theory prediction 
(see text) for Fh is given as the continuous cyan (hght) band and the rate-average of this 
prediction across the bin is indicated by the purple (dark) region. No SM prediction is 
included for the regions close to the narrow cc resonances. 



Equation ([T| is used to describe the signal angular distribution. The background 
angular and mass shapes are treated as independent in the fit. The angular distribution 
of the background is parameterised by a second-order Chebychev polynomial, which is 
observed to describe well the background away from the signal mass window (5230 < 
mK+^,+^,- < 5330MeV/c2). 

The resulting values of A^b and Fh in the bins of are indicated in Fig. [3] and in 
Table [T] The measured values of A-pB are consistent with the SM expectation of zero 
asymmetry. The 68% confidence intervals on A-pB and Fh are estimated using pseudo- 
experiments and the Feldman-Cousins technique [34|. This avoids potential biases in 
the estimate of the parameter uncertainties that come from using event weights in the 
likelihood fit or from the boundary condition (|Afb| < -Ph/2). When estimating the 
uncertainty on ApB (-^h), -^h (^fb) is treated as a nuisance parameter (along with the 
background parameters in the fit). The maximum- likelihood estimate of the nuisance 
parameters is used when generating the pseudo-experiments. The resulting confidence 
intervals ignore correlations between ApB and Fh and are not simultaneously valid at the 
68% confidence level. 



6 Systematic uncertainties 

For the differential branching fraction measurement, the largest source of systematic 
uncertainty comes from an uncertainty of ~ 4% on the -B^— > K^J/il) and J/ip — t- 



branching fractions 26 . The systematic uncertainties are largely correlated between the q 
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Table 1: Signal yield (iVsig), differential branching fraction (dS/dg^), the parameter Fh 
and dimuon forward-backward asymmetry (Afb) for the — )■ K~^n'^n~ decay in the 
bins used in the analysis. Results are also given in the Kq^ < 6GeVVc^ range where 
theoretical uncertainties are best under control. 



g2 (GeVVc^) 


N ■ 


dB/dq^ (10-8 GeV-^c^) 














0.05 - 


2.00 


159 ± 14 


2.85 ±0.27 ±0.14 


0.00 


+0.12 
-0.00 


+0.06 
-0.00 


0.00 


+0.06 
-0.05 


+0.03 
-0.01 


2.00- 


4.30 


164 ± 14 


2.49 ±0.23 ±0.10 


0.14 


+0.16 
-0.10 


+0.04 
-0.02 


0.07 


+0.08 
-0.05 


+0.02 
-0.01 


4.30- 


8.68 


327 ± 20 


2.29 ±0.16 ±0.09 


0.04 


+0.10 +0.06 
-0.04 -0.04 


-0.02 


+0.03 +0.03 
-0.05 -0.03 


10.09 - 


12.86 


211 ±17 


2.04 ±0.18 ±0.08 


0.11 


+0.20 +0.02 
-0.08 -0.01 


-0.03 


+0.07 +0.01 
-0.07 -0.01 


14.18 - 


16.00 


148 ± 13 


2.07 ±0.20 ±0.08 


0.08 


+0.28 
-0.08 


+0.02 
-0.01 


-0.01 


+0.12 
-0.06 


+0.01 
-0.01 


16.00 - 


18.00 


141 ± 13 


1.77 ±0.18 ±0.09 


0.18 


+0.22 
-0.14 


+0.01 
-0.04 


-0.09 


+0.07 
-0.09 


+0.02 
-0.01 


18.00 - 


22.00 


114 ± 13 


0.78 ±0.10 ±0.04 


0.14 


+0.31 
-0.14 


+0.01 
-0.02 


0.02 


+0.11 
-0.11 


+0.01 
-0.01 


1.00- 


6.00 


357 ±21 


2.41 ±0.17 ±0.14 


0.05 


+0.08 
-0.05 


+0.04 
-0.02 


0.02 


+0.05 
-0.03 


+0.02 
-0.01 



bins. The uncertainties coming from the corrections used to calibrate the performance of the 
simulation to match that of the data are at the level of 1 — 2%. The uncertainties on these 
corrections are limited by the size of the D*^ D^{^ K~7r^)7r~^ and J/i/j A*"*"/^" control 
samples that are used to estimate the particle identification and tracking performance in 
the data. The signal and background mass models are also explored as a source of possible 
systematic uncertainty. In the fit to the K'^fi^fi~ invariant mass it is assumed that the 
signal line-shape is the same as that of the -B"*"— )■ K^J/ip decay. In the simulation, small 
differences are seen in the mass resolution due to the different daughter kinematics 
between low and high q^. A 4% variation of the mass resolution is considered as a source 
of uncertainty and the effect on the result found to be negligible. 

For the extraction of A-pB and Fh, the largest sources of uncertainty are associated 
with the event weights that are used to correct for the detector acceptance. The event 
weights are estimated from the simulation in O.SGeV^/c"^ wide q^ bins (driven by the 
size of the simulated event sample). At low (f, the acceptance variation can be large (at 
extreme values of cos 9^) over the q^ bin size. The order of the uncertainty associated 
with this binning is estimated by varying the event weights by half the difference between 
neighbouring q^ bins and forms the dominant source of systematic uncertainty. The size 
of these effects on ApB and Fh are at the level of 0.01 — 0.03 and 0.01 — 0.05 respectively, 
and are small compared to the statistical uncertainties. Variations of the background mass 
model are found to have a negligible impact on Ape and F-^. 

The background angular model is cross-checked by fitting a template to the angular 
distribution in the upper mass sideband and fixing this shape in the fit for A-pB and Fh 
in the signal mass window. This yields consistent results in every q^ bin. Therefore, no 
systematic uncertainty is assigned to the background angular model. Two further cross 
checks have been performed. Firstly, Ap-Q has been determined by counting the number of 
forward- and backward-going events, after subtracting the background. Secondly, Fh has 
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been measured by fitting tlie |cos0^| distribution, wliicli is independent of Ap-^- Consistent 
results are found in both cases. 



7 Conclusions 



The measured values of ApB and Fh are consistent with the SM expectations of no 
forward-backward asymmetry and Fh ~ 0. The differential branching fraction of the 
— )■ K~^fi~^fj,~ decay is, however, consistently below the SM prediction at low g^. 
The results are in good agreement with, but statistically more precise than, previous 

Belle 



measurements of dB/dq and A-pB from BaBar 35 



36 and CDF 37 . Integrating 



the differential branching fraction, over the full q range, yields a total branching fraction 



of (4.36 ±0.15 ±0.18) X 10"^ 
(4.8 ±0.4) X 10"^ l26|. 



which is more precise than the current world average of 
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